Dissipative binding of atoms by non- conservative forces 
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The formation of molecules and supramolecular structures results from bonding by conservative 
forces acting among electrons and nuclei and giving rise to equilibrium configurations defined by min- 
ima of the interaction potential. Here we show that bonding can also occur by the non-conservative 
forces responsible for interaction-induced coherent population trapping. The bound state arises in a 
dissipative process and manifests itself as a stationary state at a preordained interatomic distance. 
Remarkably, such a dissipative bonding is present even when the interactions among the atoms 
are purely repulsive. The dissipative bound states can be created and studied spectroscopically in 
present-day experiments with ultracold atoms or molecules and can potentially serve for cooling 
strongly interacting quantum gases. 



PACS numbers: 32.80. Qk, 34.20.Cf, 03.65.Yz, 37.10.De 

In most experiments investigating coherent quantum 
dynamics, dissipation is an undesirable process. How- 
ever, there have been several recent theoretical proposals 
for turning controlled dissipation into a useful resource, 
e.g., for the realization of interesting many-body quan- 
tum states [IHS] . The main idea is to engineer the inter- 
action with the environment such that the combination 
of coherent and dissipative dynamics drives the system 
in question to a stationary state identical to the quan- 
tum state of interest. The feasibility of such a reservoir 
engineering has already been demonstrated experimen- 
tally [4. 

In our work, we employ a similar idea to achieve 
a dissipative bond between two interacting atoms (or 
molecules), i.e., a stationary state of the scattering dy- 
namics that confines the atoms at a fixed distance. Our 
results thus allow to extend the notion of a bonding mech- 
anism from purely conservative to dissipative forces. Ow- 
ing to the unprecedented control available for ultracold 
atoms and molecules [5 , the characteristics of the re- 
sulting molecules such as bond lengths and spectroscopic 
properties are highly tunable. 

As a specific physical implementation, we first consider 
two interacting atoms whose motion is restricted to one 
spatial dimension, see Fig. [T] (a). For ultracold atoms 
and molecules, such a constraint can be enforced by ap- 
propriate trapping potentials [5 . We assume the inter- 
nal level structure of the atoms as shown in Fig. [l] (b): 
each has one internal state |1) exhibiting a distance- 
dependent interaction shift, and two other noninter act- 
ing states 1 2) and |3), with |2) undergoing spontaneous 
decay into states |1) and |3). This particular level struc- 
ture is quite common in ultracold atoms and molecules, 
here, we focus on two different possible realizations based 
on Rydberg-dressed atoms [6HS] or laser-cooled molecules 
[9HTT]. While our approach is general, we focus on dipole- 
dipole interactions, as this brings in additional tunability 



due to the anisotropy of the interaction potential [12]. 
We generalize our results to higher spatial dimensions 
and discuss the prospects of using dissipative bonding 
as a cooling mechanism for strongly interacting quantum 
gases. 

The process behind the formation of the bond is co- 
herent population trapping (CPT) [13] [14]. In our setup, 
we include two counter-propagating laser beams; the one 
corresponding to Rabi frequency Q.- is resonant with the 
transition between states |2) and |3), and the one with 
Rabi frequency Vtj^ is detuned by A from the resonance 




FIG. 1: Setup of the system, (a) Strong confinement by 
external fields restricts the atoms' motion to one spatial di- 
mension parallel to the two counter-propagating laser beams 
represented by the driving fields with Rabi frequencies Q±; 
(b) Internal level structure of the atoms. Two internal states 
|1) and |3) are coupled to a metastable state |2) spontaneously 
decaying with a rate 7. In the non-interacting limit the field 
is detuned by A from the |1) — |3) transition; state |1) is 
subject to a dipolar interaction shift 5{r) which depends on 
the interparticle distance r; |2) and |3) are noninteracting. 
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between |1) and |3). CPT entails a dark state, i.e., a 
stationary state which cannot absorb photons. 

For a system of two atoms subject to the electromag- 
netic field treated as a reservoir [T5j, the dynamics is 
obtained from the quantum master equation. 



dp 
'di 



whose stationary solution is given by the condition 
dp/dt = 0. Here, H is the Hamiltonian describing the 
two atoms and their interaction with the laser fields, p 
is the density operator describing the two atoms, and 7^ 
and Cji are the rates and corresponding jump operators 
associated with spontaneous decay. The coherent part 
of the dynamics according to the quantum master equa- 
tion ([1]) is given by the two-atom Hamiltonian, 
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Here, z = 1,2 and k label the atoms and their corre- 
sponding momentum states, and Vdd{Q) is the Fourier 
transform of the dipole-dipole interaction potential. The 
dissipative part of Eq. ([T]) contains the rates 7n = 7 and 
jump operators = J]/, |/c + A/Cn, jn)(2, k\i^ in the Lind- 
blad form, responsible for the decay of each atom from 
state |2). The index = 1,2 runs over the two atoms, 
while jn = 1,3 accounts for the two final states, and Akn 
contains all possible values of the emitted photon's wave 
vector [16 . 

From the master equation ([T]), one can extract an ef- 
fective non-Hermitian Hamiltonian, i^eff = H — iVd^ con- 
taining a dissipative potential. 
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The zero-energy eigenstate of i^eff corresponds to the 
dark state IV^dark) if its dissipation vanishes, i.e., (c]^c^) = 
0. In the non-interacting but resonant case, 5{r) = A = 
0, the dark state for each atom is given by (|1) — |3))/a/2. 
Taking the coupling to translation via the Doppler effect 
into account retains this dark state and gives rise to the 
well-known velocity-selective CPT [14 . In contrast, here 
we consider the regime where the kinetic energy is small 
compared to the dipolar interaction between the atoms, 
which in turn is small compared to both the laser driving 
and the decay, i.e., {hk)'^/2m <C A^5{r) <C 7,^±. In the 
limit of infinite mass m, we find via perturbation the- 
ory that the dissipation has a minimum when the atoms 



are separated by a distance rd where the dipolar inter- 
action d{r) destructively interferes with the detuning A 
the most. 
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Here d is the dipole moment of state |1), Sq the vacuum 
permittivity, and h Planck's constant. While there is still 
some residual dissipation for atoms confined to the vicin- 
ity of Td^ the probability to find the atoms separated by 
Td is strongly enhanced. Such a confinement of the atoms 
amounts to the formation of a dissipative bond. Addi- 
tional perturbations arising due to a finite mass lead to a 
stationary state that exhibits a distribution of distances 
sharply centered around Td rather than a single fixed dis- 
tance. 

We now turn to numerical simulations to confirm the 
existence of the dissipative bond. We solve the quan- 
tum master equation, Eq. ([T]), using the wave- function 
Monte-Carlo method [16]. During the time evolution the 
atomic positions are redistributed by the kinetic energy 
and photon recoil, and the atomic population is accu- 
mulated in the vicinity of Td on a fast timescale of a few 
tens of 7~^. For an ensemble of atoms with randomly dis- 
tributed initial conditions we observe the appearance of a 
quasi-stationary distribution strongly peaked around r^/, 
see Fig. [2] The results presented in Fig. [2] correspond to 
a state after the evolution time of 75 7"^, averaged over 
250 realizations with random initial conditions. In each 
realization two atoms are confined in a one-dimensional 
box L ~ 6rd^ which corresponds to the three-dimensional 
particle density of 3 x 10^^ cm~^. 

For even longer times the spatial distribution of the 
atoms appears to evolve towards a completely delocalized 
equal- weight superposition of positions compatible with a 
relative distance rd ^ 500 nm. This constitutes evidence 
that in the long-time limit the state corresponding to the 
dissipative bond is a nearly pure state with a low entropy. 
Once the dissipative bond has been formed, it is possible 
to further reduce the probability of subsequent sponta- 
neous decay events by decreasing the Rabi frequencies 
Q±. This leads to an effective decay rate 7eff = 7^2, 
with p2 the total probability of any of the atoms being in 
state 1 2). We find that the reduction of the Rabi frequen- 
cies does not lead to a significant increase of the spread 
about rd^ and that lifetimes of the bound state exceeding 
7eff^ ~ 0.1 s can be achieved. 

Experimental observation of the resulting bound states 
can be performed by making use of the techniques well- 
established in the area of ultracold quantum gases. In 
particular, a suitable setting would be a many-body sys- 
tem, in which the atoms are essentially uncorrelated 
at the beginning, with an average interatomic distance 
much larger than r^^, so that many-body effects could be 
neglected. Then, the existence of the bound state will 
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^ (j^j^^ Relative distance (nm) 

FIG. 2: Probability density distributions indicating the formation of the dissipatively bound states: (a) the absolute coordinates 
of the two particles, ri, r2, are delocalized; (b) the projection onto the relative distance r = |ri — r2| exhibits a sharp peak 
around rd = 500 nm. Parameters correspond to a pair of Rydberg-dressed Cesium atoms (see text). 



appear as a sharp peak in the pair-correlation function, 
see Fig. [2] (b) , which is readily accessible through Bragg 
scattering [17 or noise correlation spectroscopy p!8l [T9] . 

The above results can be generalized to higher spatial 
dimensions, thereby allowing to study the motion of the 
bound atoms. This merely requires the inclusion of ad- 
ditional counter-propagating laser beams as used in con- 
ventional laser cooling setups [20 . For a two-dimensional 
configuration, we consider the case when the quantization 
axis defined by the electric dipole moment is perpendic- 
ular to the plane of the motion, i.e., when the dipolar in- 
teraction is always repulsive. Then, the dissipative bond 
restricts the relative radial motion to a ring where the 
atoms are separated by Vd. The effective Hamiltonian 
i^eff for the remaining angular motion is then described 
by a two-dimensional rigid rotor, i^eff = BJ'^^ with 
the projection of the angular momentum on the quanti- 
zation axis. The rotational constant B = /(mr'^ can 
be widely tuned by changing Vd] for typical experimen- 
tal realizations, we find that values oi B ^ 100 — 1000 Hz 
are realistic, i.e., by far exceeding the effective decay rate 
7eff. The rotational spectrum of the bound atoms can be 
then probed with standard techniques of molecular spec- 
troscopy. 

In three spatial dimensions, the constrained relative 
motion of the two atoms is even richer due to the 
anisotropy of the dipolar interaction. The dynamics is 
confined to a surface defined by 

=r^(l-3cos^i?), (5) 

with the angle between the collinear dipoles and the 
interatomic radius-vector. Due to the rotational sym- 



metry about the quantization axis, the projection of the 
angular momentum, Jz^ is conserved, allowing to reduce 
the problem to an effective one-dimensional problem. 
The effective potential corresponds to an azimuthally- 
symmetric double well centered at ^ = 7r/2, resulting in 
an almost equidistant spectrum consisting of tunneling 
doublets. Fig. [3] (a), similar to those encountered in mi- 
crowave spectra of ammonia [21] or molecules in strong 
laser fields [22]. Fig.|3](b) features the probability density 
distributions on the three-dimensional surface. Although 
the results shown in Fig. [3] correspond to = 0, energy 
spectra and wavefunctions for other values of Jz exhibit 
similar features. 

Experimental implementation. — We consider two dif- 
ferent realizations of our setup depicted in Fig. [l] 
based on Rydberg-dressed Cesium atoms and laser- 
cooled molecules. In the case of Rydberg dressing |6j- 
8 , the states |1) and |3) are chosen as different hyperfine 
components of the electronic ground state, 6^ 51/2 , of Ce- 
sium. The state |1) is provided with a dipole moment of 
(i = 15D perpendicular to the ID trap due to Rydberg 
dressing in an external electric field; laser fields Vt± =7/2 
drive the 6^5'i/2 ^ ^^P?,/2 transition having a linewidth 
of 7 = 27r X 5.2 MHz. We would like to stress that in con- 
trast to the previous Rydberg dressing proposals for the 
observation of interactions in a Bose-Einstein condensate, 
the requirements on coherence times for the observation 
of the dissipative bond are much less stringent; preserv- 
ing coherence on the timescale of the rotational constant 
B 100 Hz can be achieved in present-day experiments 
[23H28] , see Ref. [29] for a recent review. 

Although we exemplified our scheme by the case of 
bonding between a pair of atoms, it is equally applicable 



4 




FIG. 3: Three-dimensional eigenstates [Jz — 0): (a) lowest energy levels featuring tunneling doublets due to the azimuthal 
symmetry of the interaction potential; (b) probability density distributions. The angle tl) gives the orientation of the collinear 
dipoles with respect to the vector connecting them. The quantum number n gives the number of nodes of the probability 
density distribution (purple color) depending on the angle i^. 



to bonding between a pair of molecules, such as strontium 
monofluoride (SrF) that is being laser-cooled ^Oj. For 
SrF, different hyperfine components, F = 0,1, of the 
X^i;+('U = 0; = 1, J = 1/2) state are chosen as |1) and 
1 3). Here, a dipole moment of 3.5 D in the rotating frame 
is imprinted on |1) via microwave dressing [30 , while the 
fields Vt± = 7/2 drive transitions to the electronically 
excited A'^Ui/2{v' = 0;N = 0,J' = 1/2) whose natural 
linewidth is 7 = 27r x 7 MHz. 

Finally, we note that the mechanism underlying the 
formation of the dissipative bond could be used for cool- 
ing of strongly-interacting many-body systems. Clearly, 
the appearance of a low entropy stationary state inde- 
pendent of the initial conditions already amounts to a 
demonstration of strong cooling. In the high-density 
regime, the correlation length rd will be modified by an 
effective coordination number, embodying the interac- 
tions of an individual atom with its surroundings in the 
stationary state. Despite this renormalization, the for- 
mation of dissipatively bound complexes [31^ in free space 
or even the realization of dipolar crystals |32l [33] appear 
possible. 
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